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Abstract. We introduce a theoretical model to scrutinize the conductivity of small 
polarons in one-dimensional disordered systems, focusing on two crucial -as will be 
demonstrated- factors: the density of states and the spatial extent of the electronic 
wave function. The investigation is performed for any temperature up to 300 K and 
under electric field of arbitrary strength up to the polaron dissociation limit. To 
accomplish this task we combine analytical work with numerical calculations. 
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1. Introduction 

The density of states (DOS) is the heart of any physical system in the sense that 
its structure and magnitude crucially affect all physical properties. Particularly, this 
holds for the response of charge carriers to external stimuli such as electric or magnetic 
fields, temperature gradients or temperature variations, i.e. the transport properties. In 
disordered materials, the random distribution of their constituents drastically affects the 
character of the carriers and the transport mechanisms. Under certain circumstances, 
the presence of disorder induces carrier localization and hopping becomes the chief 
transport mechanism. Hence, the electronic wave function spatial extent (a -1 ), being a 
measure of the carrier localization, becomes a parameter of vital importance. 

One-dimensional (ID) systems have been recently considered to be among the 
most promising materials for nanotechnology. In particular, an increasing amount of 
experimental and theoretical work has been devoted to the electrical properties of ID 
amorphous semiconductors, amorphous carbon, doped polymers, conjugated polymers 
and organic materials [Il|2i|3ilIJ[5j[6j[7J|8j[9j [HJJ [EE] , [121 and references therein] . 

Given that DNA has been placed among the most promising organic materials 
for nanotechnology, Triberis et al. [13], studied DNA as a ID disordered molecular 
"wire" in which small polarons are the charge carriers. Based on the Generalized 
Molecular Crystal Model (GMCM) [2] and theoretical percolation arguments, they 
studied small polaron hopping along the DNA double helix and in the presence of 
low electric field (F). Ignoring the effect of correlations, an analytical expression for 
the strong temperature (T) dependence of the electrical conductivity (cr) was obtained 
which reproduced the experimental data reported for A-DNA [15J and for poly(dA)- 
poly(dT) DNA [16] at high temperatures. The theoretical analysis also permitted the 
evaluation of the maximum hopping distance and its T-dependence, supporting the idea 
of multi-phonon assisted hopping of small polarons between next nearest neighbors of 
the DNA molecular "wire". Taking into account the effect of correlations (cr), Triberis 
and Dimakogianni (T2J [17] showed that lna cr oc T _1//2 holds for high as well as for 
low temperatures. This reproduced the strong o~(T) at high temperatures reported for 
A-DNA [151 H] an d poly(dA)-poly(dT) DNA [16], while, including correlations, the 
evaluation of the maximum hopping distance led to systematically longer values than 
those evaluated ignoring correlations [13], supporting experimental evidence for long 
range charge migration along the DNA double helix [T9], [20], |2~Tj . 

In addition, even under moderate electric fields, strong nonlinearities of c(F) in 
ID disordered systems have been observed. In the variable range hopping regime and 
at low temperatures, Fogler and Kelley [22] investigated theoretically the effect of a 
finite electric field on the resistivity. They took into account the existence of highly 
resistive segments (breaks) on the conducting path of the carriers in ID systems and 
found that the role of the breaks diminishes and eventually becomes insignificant as F 
increases. Ma et al. [23] described hopping transport and the conductivity of ID systems 
with off-diagonal disorder. Investigating the T-dependence of the hopping conductivity, 



Density of states and extent of wave function 



3 



they showed that it increases with the increase of T taking much larger values than in 
the case of the Anderson model with pure diagonal disorder. They also studied the in- 
dependence of the conductivity to find that at low F the hopping conductivity conforms 
with the ohmic law, but at strong fields it presents non-ohmic characteristics. 

Triberis and Dimakogianni [21] studied the cr(T, F) behaviour under the influence 
of moderate electric fields up to ~ 10 5 Vm" 1 , when small polarons are transported in a 
disordered ID environment, at high and low temperatures. The analytical expressions 
obtained for a(F, T), were applied to experimental findings concerning charge transport 
in polydiacetylene quasi-lD single crystals [9]. It was shown that at low electric fields 
the hopping conductivity conforms with the ohmic law while increasing the electric 
field the conductivity presents non-ohmic characteristics. The transition from the 
ohmic to the non-ohmic behaviour starts for smaller values of F at lower temperatures 
and the rate of the increase of a is greater the lower T is. These conclusions were 
in a qualitative agreement with theoretical results referred to variable range hopping 
[221 1231 125] . Dimakogianni and Triberis [2E] also investigated the effect of correlations 
on the non-ohmic behaviour of the small polaron hopping conductivity in ID and 3D 
disordered systems. They concluded that the inclusion of correlations results to a much 
stronger dependence of the conductivity on the magnitude of the applied electric field 
compared to the uncorrelated case. The deviation of the conductivity from the ohmic 
behaviour appears twice as fast when correlation effects are taken into account, for a 
given applied electric field as the temperature increases. 

In the present work, taking into account the directionality imposed by the electric 
field on the transport path of the carriers, we examine the role of the magnitude of the 
density of states and the extent of the electronic wave function and calculate a. The 
aim of the present work, is to investigate o~(T, F) for all reasonable T and F values, i.e. 
from 10 up to 300 K and up to the F values where polarons cease to exist. This is done 
varying the density of states by orders of magnitude around values which are relevant 
to common ID systems [271 1211 [29j [30] and varying the extent of the electronic wave 
function from 1 to 5 A, i.e. reasonable values for common organic molecules [3X1132"] . We 
demonstrate that F plays both a constructive energetic role by offering energy for the 
carrier jumps and simultaneously a destructive role, in the sense that the stronger it is 
the more it forces the polaron to jump opposite to the F direction prohibiting forward 
jumps to neighboring sites. 

In Section [2] we present our theoretical model including the basic analytical 
expressions at high and low temperatures. According to the mathematical analysis of 
the Generalized Molecular Crystal Model [HJ [33], it is the condition huo <C ksT 
(hujQ S> ksT) that determines the high (low) temperature regime. This mathematical 
analysis leads to the evaluation of the intrinsic transition rate, which differs at high 
temperatures (multi-phonon assisted hopping), compared to that at low temperatures 
(few-phonon assisted hopping). Which temperature range in real systems is indeed high 
or low depends on the system under study. Our numerical results together with the 
relevant discussion are staged in Section [3j In this way we examine the conductivity of 
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small polarons in one- dimensional disordered systems, and demonstrate that the density 
of states and the spatial extent of the electronic wave function are two crucial factors 
for its behaviour. The temperature and the electric field ranges that we consider are 
very broad. In particular, the electric field is varied from very low up to the polaron 
dissociation limit (~ 1 x 10 8 Vm" 1 ). Finally, in Section H] we state our conclusions. 

2. Theory 

2.1. Generalized Molecular Crystal Model 

In the context of GMCM we consider a ID deformable "wire" consisting of "molecular 
lattice sites" across which small polarons are transported in the presence of disorder. 
By €i(0), and e,-(0) we denote the energies of an electron on site at vector positions 
and Tj, respectively, if the "molecular lattice sites" are constrained not to be displaced 
in response to the presence of the electron. Due to the disorder these local electronic 
energies, ej(0), and 6^(0) are not equal. The energetic non-equivalence of the two sites 
will affect the small polaron's binding energy, E^(i), in the sense that, the lower the 
local electronic energy is the more localized the electronic wave function will tend to be 
and consequently the larger its binding energy will be. Assuming that the stiffness of 
the "molecular lattice" is unaltered, the difference in binding energy means a difference 
in the electron- lattice interaction parameters A» and Aj i.e. -Ej(xj) = 6j(0) — AjXj and 
Ej(x.j) = ej(0) — A/Xj with Ai ^ Aj. Here, E%{xi) is the electronic energy of the 
system of the electron and the isolated molecule with configurational coordinate x«, 
which represents the deviation of the atoms of the molecule at position rj from their 
equilibrium configuration i.e. the local vibrational displacement coordinate. 

The GMCM p3] is based on a generalized "hopping model" Hamiltonian of the 

form 

< m\H\n >=< m\H + V\n >= £i> k }(%5{n k },{n k ,}+ < m\V\n >, (1) 

the < m| V|ra > term [14] is the overlap part of the Hamiltonian, |n >= \i, {n^} > are 
the eigenstates of H, and Hq is the zeroth-order (i.e. for electronic overlap integral of 
the tight-binding theory J=0) Hamiltonian with corresponding eigenvalues 

#i,{n k } = £i(0) - E b (i) + £ hu k {n k + \). (2) 

k 1 

Here, {n^} represents the totality of the vibrational quantum numbers (...., rik, ■••) for 
the occupation of the site with position vector r^, and 

W) = ^E(4 2 /2M^), (3) 
iV k 

is the small polaron binding energy. iV is the number of "molecular lattice sites" and 
M is the appropriate reduced atomic mass. The relation between and its associated 
wavevector k, i.e. the dispersion relation, is given by: 

u£ = u% + ulY, coa ( k ' h! )> ( 4 ) 

k 
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where |k| = 27rp/N, the integer p lying in the range — (iV — 1)/2 < p < (N — l)/2, and h' 
indexes the nearest neighbors (r-j+h') of an arbitrary site r$. uq is the harmonic oscillator 
frequency associated with the configurational coordinate of the isolated molecule. The 
relation u\ <C Uq determines the weak dispersion limit. 

Equations (2) and (3) show the essential features of the GMCM which are: 

1. site-dependent local electronic energy ej(0). 

2. site-dependent electron-lattice interaction parameter, and concomitant binding 
energy E b (i). 

The knowledge of < m|V^|n >, permits the evaluation of the "microscopic" small 
polaron velocity operator [Ml 135] , 

iijj =< m|u|n >= (— ) < m|V|n > (r,- — r^), (5) 

the charge current density operator, 

jij = n c quij, (6) 

where n c is the charge carrier concentration, and q is the carrier's charge, and thus the 
"microscopic" electrical conductivity [36J, 

roc r/3 

CTij= alt dp < j(-ihp)j(t) >, (7) 
Jo Jo 

where (3 = l/fc^T. The mobility, /j,^, and consequently the diffusion constant, given by 
Dij = fiij/e(3, are determined and lead to the "microscopic" jump rate which reads: 

D tl _ 

12 



^ = 7^772 (8) 



r ii 



Assuming that the dependence on the spatial separation Rij, of the two sites is 
exp(— 2aRij), the "microscopic" intrinsic transition rate, 7^, for a small polaron hopping 
from a site % to an empty site j is given by 

jij = exp(-2aRij)Lij. (9) 

The treatment refers to the non(anti)-adiabatic limit, i.e. in the physical situation 
where the electron is no longer able to follow rapid fluctuations of the lattice and, hence, 
it does not respond quickly enough to the occurrence of a coincident event in order to 
overcome the energy barrier. In this case, J can be treated as a small perturbation in 
the lowest order [Ml EHl [39]- 

The expansion of the model to include the influence of possible strong local 
interparticle correlations might be interesting as intercarrier interactions exist in real 
systems. However, this is beyond the aim of the present work. 

2.2. Hopping at high temperatures 

When a carrier hops from site i of energy Ei to site j of energy Ej, at a distance R^, 
the intrinsic transition rate between the two localized states at high (h) temperatures 
(hu <C ksT [H]) is 

7S = 75 expt-2^) exp (- jSL) X { Jg>. I > (10) 
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Figure 1. High temperatures. (I) Left panel. Contour of constant 5ft' 1 for an initial site 
of U'. All possible final sites for the carrier, A/'(5ft' 1 ), lie on or within the contour 5ft' 1 , 
for a particular 9. (II) Right panel. Contour of integration 5ft' 1 = 5ft„„, for an initial 
site of E\ and particular 6, for the evaluation of R'p. 



Here a -1 is the spatial extent of the electronic wave function, 7q = ( J 2 /^)(7t/4£ 2 A;bT) 1 / 2 
and e 2 = [E b (i) + E b (j)]/A. E b {i) and E b (j) is the small polaron binding energy for sites 
i and j, respectively. Hence, 7^ (as well as 7^) have both spatial and energy dependence 
[HI [33]. The spatial dimensions of the system and the number of energies involved in 
the expression of the intrinsic transition rate can be considered as the coordinates of a 
"hopping space" in which the small polaron transport occurs under the influence of F. 
In this "hopping space", the most probable hop for a carrier on a site at energy E^ is 
to the empty site at closest range, i.e. to its nearest neighbor site. The average nearest 
neighbor range in the "hopping space", R nn , determines the conductivity of the system 
[40J. Thus, to evaluate the electrical conductivity we have to calculate this quantity 
first. Then, taking into account that in real space greater real forward distances will be 
hopped in the downfield direction rather than upfield, an average real forward distance 
hopped should be evaluated which, as will be presented in the following (cf. Eqs. [33ll34|) . 
leads to the mobility of the carriers and finally the overall conductivity of the system. 

From the expression of the intrinsic transition rate between two sites i and j, we 
define the range 9f& between the sites in the "hopping space" 



Taking the energies of the carrier to be mainly polaronic [14] , and using for convenience 
the terms Et and Ej instead of E b {i) and E b (j) respectively, we obtain e 2 = (Ei + Ej) /4. 
Therefore 




(11) 



K?, = 2aR ij + 



3Ej — E{ 
4k R T 



(12) 
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^ - ~^ ^ - ^ + X' cos 0, (14) 



[/' = = -e; + cos e, (is) 



where £7j > Ei for absorption and Ei/3 < Ej < Ei for emission of phonons. Introducing 
the dimensionless coordinates i2L = 2aRij, E[ = Ei/2k B T and £j = Ej/2k B T, 

where Ej > E[ for absorption and E[j3 < Ej < E\ for emission. 

Under the influence of an externally applied electric field the actual energy of the 
hop is modified [ID] 

2 3 2 1 2 3 2 1 2 
where = eF/2aksT and 9 is the angle between the directions of R\j and F. Defining 
the reduced initial ([/') and final (£") coordinates in the "hopping space" 

-E' E' = -E' + £ 
2 4 2 J 2 

the range between two sites in the "hopping space" becomes 

^ h = R! + E' - U', (16) 

where E' > 3U' for absorption and U' < E' < 3U' for emission. The indices from 9ft^ 
and R[j have been dropped. 

For the evaluation of the average nearest neighbor range, 3ft^ n , firstly we have to 
evaluate the number of unoccupied sites within a range 3l h of a particular site of U' 
as a function of T and F, Af(3l h ). The three-dimensional "hopping space" can be 
represented, for a particular 9 by a two-dimensional diagram (Fig. [Tj (I)). 

For hops of range less or equal to ^t h from an initial site of U', the final sites 
will lie on or within the contour dt h , for a particular 9, i.e. in the space defined by 
U' <E' < {n h + U')-R' and < R' < M h (Fig. Q] (I)). Thus, using Eqs.QH the number 
of empty sites enclosed by the contour Sft' 1 is 

X(K h ) = W / iV(4)[l - fiE'^^-dE'jdR'. (17) 
" Jo ./££ J J a J 

CC = \{E[ - /3R' cos 9), UC = |[2K h + E[ - 2R'(1 + § cos 9)} and 9 = 0, vr as we have 
chosen to put F exactly along the ID axis where transport takes place. N(E'j) is the 
density of states, and f(E'j) the Fermi-Dirac distribution. We take the Fermi energy 
Ep = 0. Assuming a constant density of states, N{E'-) = N(E' i ) = N , 

/•5R' 1 1 + M2-R h +E'-2R\l+t cose)) 

Arm h ) = cy in—- — 2 , p , — aN — da', (is) 



0.7T 



where C = 

2a 
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Figure 2. Low temperatures. (I) Left panel. Contour of constant 3fJ' for an initial site 
of U' . All possible final sites for the carrier, A/"(5R Z ), lie on or within the contour 5ft', 
for a particular 9. (II) Right panel. Contour of integration 5ft' = 5f nn , for an initial 
site of E* and particular 0, for the evaluation of R' F . 



2.3. Hopping at low temperatures 

The intrinsic transition rate at low (I) temperatures {Tiujq 3> k B T [53]) is 



7L = Toexp(-2ai? i ,) { f^S^ ^ > E J . (19) 



1, £j < E; 



TV J ex f — \ 

= T! ~ -] 2 [(4g2/^o) Aij/ ^°/( A ij/^o)!] ■ Following the same methodology 

as for the high temperatures, we assign every hop of the carrier to a hop in a three- 
dimensional "hopping space" defined by one spatial and two energy coordinates. 

From the expression of the intrinsic transition rate between two sites i and j, we 
define the range 3?^- between the sites in the "hopping space" 

= { If-' + E J > * . (20) 

Introducing the dimensionless coordinates R'^ = 2aRij, E* = Ei/k B T and E* = 
Ej/k B T, 

!t?' _L Z7 1 * Z?* Z?* \ Z7 1 * 

R' ip e*<e; ■ [ZL) 

Under the influence of an externally applied electric field the actual energy of the 
hop is modified [10] 

Ej - E* ~^E*- E* + m i3 cos 9, (22) 

where (3 = eF/2ak B T. Thus, 

^ (R'(l + /3cose)+E*-E*, E*>E;-PR'cos6 

\R', E* <E*-/3R' cos 9 • 1 ' 
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The indices from 3?' ■ and R'^ have been dropped. Defining the reduced initial (IF) and 
final (E') coordinates in the "hopping space" 

U* = E*, E* = E* + 0R' cos 9, (24) 

the range between two sites in the "hopping space" becomes 

(25) 

For the evaluation of the average nearest neighbor range in the "hopping space", 9ft ran , 
firstly we have to evaluate the number of unoccupied sites within a range ^R 1 of a 
particular site of U* as a function of T and F, Af(JR l ). The "hopping space" can be 
represented, for a particular 9 by a two-dimensional diagram (Fig. EJ^I)). 

For hops of range less or equal to 3t l from an initial site of U* , the final sites 
will lie on or within the contour 9^, for a particular 9, i.e. in the space defined by 
-oo < E* < (R l + U*) - R' and < R' < $t l (Fig. EI (I)). Thus, using Eqs. HQ the 
number of empty sites enclosed by the contour 3t l is 

= E / / N(E*)[1 - f(E*)]^-dE*dR'. (26) 
^7 Jo Jcc J J 2a J 

CC = — oo, UC = + E* — R'(l + (3cos9). Assuming a constant density of states, 
N(E*) = N(E*) = N , 

jV($)=CY i J ln[l+exp(3f^ / + J E;-i? , (l + /3cos0))]d J R , , (27) 



0.7T 



where again C = 



2.4- Conductivity 

We define E" = E[ for high temperatures or E" = E* for low temperatures. The 
knowledge of the number of unoccupied sites within a range %t h ' 1 for either high or low 
(h/l) temperatures, J\f(JR h ^ 1 ), permits the evaluation of the average nearest neighbor 
range, , when the carrier resides on a particular site of E", as a function of T and 
F @D] 

t!l = f ® h/l ^§T- eM-m h/l m h/l , (28) 
or equivalently 

M h Jl = [°° exp\-N($t h / l )]d$l h/l . (29) 
Jo 

The evaluation of 9ft^J, gives the range in the three-dimensional "hopping space" where 
a nearest neighbor exists that can host the carrier when the carrier hops from an initial 
site of E". However, it gives no information on the direction of the hop of the carrier. 

Considering, all sites of initial E" and assuming that all hops from these sites are 
all hops of range K^J, then in real space, these hops will be in random directions, but 
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for a hop to final sites of the same energy, greater real forward distance will be hopped 
in the downfield direction rather than upheld. Thus, summing over all final sites, for 
initial sites of E", there will be associated an average real forward distance hopped [40] 



Rf = (30) 
lot 

For high temperatures, the distance R F is evaluated by averaging R' cos 9 over the 
contour 3ft nn = const. (Fig. [TJ (II)), and hence 

R F = (31) 

The integrals I± and I2 are given in the Appendix. 



For low temperatures, the distance R F is evaluated by averaging R' cos 9 over the 
contour 9ft nn = const. (Fig. [2] (II)), and hence 

Tf F = £±£. (32) 

+ -<4 

The integrals ii, I2, I3, I a are given in the Appendix. 

Either for high or for low temperatures, having calculated the distance R' F l and 
considering that the probability of all hops is exp(— 3?^), the average rate of transport 
of carriers is v v frR! h F l exp(— Here, 1^ is a hopping attack frequency of the order 
of a phonon frequency, assumed the same for all hops. 

The mobility for small polarons of E" = E[ for high temperatures or of E" = E* 
for low temperatures reads 

~~Dl h / 1 

K4') = -f^eM-t!n), (33) 
and the conductivity of the system is 

a h '\F, T) = - eN(E")f{E")n(E")k B TdE". (34) 

j —00 



3. Results and discussion 



In the following, based on the theoretical analysis presented above, we calculate 
numerically the electrical conductivity varying the density of states and the spatial 
extent of the localized electronic wave function. Our numerical results refer for simplicity 
reasons to a constant density of states, although, typically, in ID systems the density of 
states has a strong energy dependence. One could alternatively use an energy dependent 
model for the density of states [3U H2] which is expected to influence somehow the 
conductivity. This is beyond the scope of the present paper, but could be numerically 
examined in the future via the same approach, as it is evident from Eq. [341 

We consider the range T = (160 - 300) K as high temperatures and the range 
T = (10 - 150) K as low temperatures. We investigate the influence of an electric field in 
the range F = (5 x 10 3 - 1 x 10 8 ) Vm" 1 . We did not consider higher values of F, because 
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it is generally accepted that the value of the highest electric field that the polaron can 
sustain is about 1 x 10 8 Vm" 1 @31 EHJ EES]. We take v ph = 10 12 s _1 . 

We vary the density of states by orders of magnitude around values which are 
relevant to common ID systems [27J [281 1221 [30] and the extent of the electronic wave 
function from 1 to 5 A, i.e. for sensible values for common organic molecules (3T| 132] . 
Specifically, in Figs. El HI [5] and Figs. El [91 [10]we keep a -1 = 2 A, while in Figs. El [7] we 
vary a~ l in the range 1 - 5 A. However, to keep our results as general as possible, we 
do not make any reference to a specific material. 



3.1. High temperatures 

For the temperature range T = (160 - 300) K, Fig. [3] presents a as a function of T -1 / 2 
for different magnitudes of the DOS, i.e. for N = (10 9 - 10 13 ) eV _1 m _1 [cases (a) to 
(e), respectively]. a~ l = 2 A. F = (5 x 10 3 - 1 x 10 8 ) Vm -1 . We depict our results as a 
function of T -1 / 2 to compare them with the analytically obtained formula In a oc T -1 / 2 
which holds for low up to moderate electric fields and was previously obtained by two 
of us [26] following a different theoretical treatment and taking into account the effect 
of correlations, namely 



a h ' cr (F } T) cxexp 



h.crx 1 // 2 F 2 \V2' 



(35) 



V T J V 9(T), 

Here g(T) = (2aksT / e) 2 and Tq' ct = lASa/ksNo. We observe that higher density of 
states leads to higher conductivity in such a way that is smaller for higher densities 
of states. As shown in Fig. El in this specific example, augmenting DOS by four orders 
of magnitude, the conductivity rises by approximately eleven orders of magnitude. 

From these results we realize that the T~ 1//2 -behaviour of a [26] holds for low up 
to moderate electric fields. For the smaller densities of states, i.e. for iVo = 10 9 and 
10 10 eV~ 1 m~ 1 , 4p; < i.e. the conductivity is larger for lower electric field strengths. 
This is due to the competitive role of the directionality imposed by the electric field and 
the temperature. This directionality affects destructively a when not many sites are 
available for the carrier i.e. for small densities of states. Here we notice that according 
to Eq. [16] and Fig. CEp), the electric field affects the range between two sites in the 
"hopping space" both for the absorption and the emission branch. On the contrary, 
this effect does not appear at low temperatures (T = (10 - 150) K), discussed later on 
Subsection 13.21 (Fig. [8]) because in the corresponding expression for the range between 
two sites in the "hopping space" at low temperatures (Eq. [251 and Fig. 12(1)), the electric 
field affects only the finite area of the absorption branch. We mention that the electric 
field plays a constructive role, too, due to its energy offer to the carriers. At N = 10 11 
eV _1 m _1 it seems that the available sites are numerous enough so that the directionality 
of the electric field hardly affects the conductivity. For higher densities of states, i.e. 
for iVo = 10 12 and 10 13 eV _1 m _1 , only the constructive energetic influence of the 
electric field appears. Now ^ > i.e. the conductivity is larger for higher electric 
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Figure 3. a versus T 1 / 2 , varying the density of states i.e. for No — 10 9 - 10 13 
cV-W 1 [cases (a) to (e), respectively]. oT x = 2 A. T = (160 - 300) K and 
F = (5 x 10 3 - 1 x 10 8 ) Vm" 1 . 



field strengths. Finally, we observe that the temperature has a greater effect on the 
conductivity, the smaller the density of states is. Another aspect of the behaviour of 
the conductivity for different DOS is shown in Figure HJ iVo = 10 9 - 10 13 eV _1 m _1 , 
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Figure 4. a versus F, varying the density of states i.e. for Nq = 10 9 - 10 13 eV 1 m 1 
[cases (a) to (e), respectively], or 1 = 2 A. T = (160 - 300) K and F = (5 x 10 3 - 
1 x 10 8 ) Vm" 1 . 



T = (160 - 300) K and F = (5 x 10 3 - 1 x 10 8 ) Vm- 1 . 

Let us denote by cr the ohmic value of the conductivity, i.e. a = lim^_ ) .o( cr )- 
In order to show the deviation of o from <7o under the influence of both F and T we 
present Fig. |5]which shows a/ao versus (3 = eF/2akBT for (a) N = 10 9 eV _1 m _1 and 
(c) iVo = 10 11 eV^m- 1 , respectively. F = (5 x 10 3 - 1 x 10 8 ) Vm- 1 and T = (160 - 300) 
K. We observe that the effect of the external stimuli F and T on a depends strongly 
on the value of the density of states that characterizes the system. For (a) iV = 10 9 
eV _1 m _1 the conductivity decreases from its ohmic value in the specific range of F and 
T, while for (c) iVo = 10 11 eV~ 1 m~ 1 the conductivity generally increases and it is higher 
for higher temperatures. We notice that the variation of a /o~q versus (3 = eF/2akBT is 
generally small especially in contrast to the corresponding variation at low temperatures 
studied later on Subsection 13.21 

Figure |6] shows the conductivity for different values of the spatial extent of the 
localized electronic wave function, i.e. for a -1 = (1 - 5) A. Here we have chosen case 
(c) iVo = 10 11 eV^m" 1 for the density of states. T = (160 - 300) K and F = (5 x 10 3 - 
1 x 10 8 ) Vm" 1 . We observe that smaller a -1 (more localized carriers) leads to smaller 
a. Particularly, five times increase of a~ l leads to two orders of magnitude greater 
conductivity. In addition, in Fig. [7] we observe that 4| > for any temperature when 
a' 1 = 3, 4 and 5 A, while g < when a' 1 = 1 A. For the case a^ 1 = 2 A cf. Fig. Efll). 
In other words, the strength of the localization which determines the size of the formed 
polaron, along with the density of states which characterizes the system, are both two 
key factors for the conductivity and its dependance on F and T. 
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Figure 5. a/a versus /3 = eF/2ak B T. (I) Left Panel. Case (a) N a = 10 9 eV -1 m -1 . 
(II) Right Panel. Case (c) A^ = 10 11 eV -1 m _:t . In both panels aT 1 = 2 A, F = (5x 10 3 
- 1 x 10 8 ) Vm- 1 and T = (160 - 300) K. 
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Figure 6. (I) Left Panel, a versus T- 1 ' 2 . (II) Right Panel, cr versus F. The spatial 
extent of the localized electronic wave function is taken aT 1 = (1 - 5) °A. Here we have 
chosen case (c) Nq — 10 11 eV~ 1 m _1 for the density of states. T = (160 - 300) K and 
F = (5 x 10 3 - 1 x 10 s ) Vm- 1 . 



3.2. Low Temperatures 

For the temperature range T = (10 - 150) K, Fig. [8] presents a as a function of 
for different magnitudes of the DOS, i.e. iVo = (10 9 - 10 11 ) eV _1 m _1 [cases (a) to (c), 
respectively], a" 1 = 2 A. F = (5 x 10 3 - 1 x 10 s ) Vm^ 1 . Again, we depict our results 
as a function of T _1//2 in order to compare them with the analytically obtained formula 
lncr oc which holds for low up to moderate electric fields and was previously 

obtained by two of us [26] following a different theoretical treatment and taking into 
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Figure 7. a/ao versus j3 — eF/2aksT for a 1 = (1 - 5) A. For the density of states 
we have chosen (c) N = 10 11 eV -1 m _1 . T = (160 - 300) K and F = (5 x 10 3 - 1 x 10 8 ) 
Vm" 1 . 



account the effect of correlations, namely 



a l > cr (F,T) ocexp 



1/2' 



T J \ g(T) : 



(36) 



Here Tq CT = 1.96a/ ks^Q. We observe that higher density of states leads to higher 
conductivity in such a way that is smaller for higher densities of states. In total, 
in this specific example, augmenting DOS by two orders of magnitude increases the 
conductivity by approximately tens of orders of magnitude. From these results we 
realize that the T~ 1//2 -behaviour of a [26J holds for low up to moderate electric fields. 
For higher values of F the conductivity deviates from the T _1 / 2 -behaviour as T decreases 
and this deviation appears to be larger the stronger the electric field is. This deviation 
is due to the constructive energetic contribution of the electric field which leads to 
the increase of the number of available sites that can host the carrier, i.e. the range 
9ft nn does not depend solely on T. As the temperature further decreases, for strong 
enough electric fields, the range 5i nn depends exclusively on the applied electric field, 
as essentially all hops are downward in energy. As a result, the conductivity does not 
depend on the temperature. In other words, there is a transition from thermally-assisted 
to field-assisted hopping. 

We remind the reader that in the high temperature range T = (160 - 300) K 
discussed in subsection 13.11 the electric field affects the range between two sites in the 
"hopping space" both for the absorption and the emission branch, leading also to the 
appearance of the destructive role of the electric field. In contrast, this effect does not 
appear here in the low temperature range T = (10 - 150) K, because in the corresponding 
expression for the range between two sites in the "hopping space" the electric field affects 
only the finite area of the absorption branch. 
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Figure 8. a versus T" 1 / 2 for cases (a) N = 10 9 eV -1 m - \ (b) N = 10 10 eV^m" 1 
and (c) iV = 10 11 eV^m -1 . a^ 1 =2 A. T = (10 - 150) K and F = (5 x 10 3 - 1 x 10 8 ) 
Vm _1 . 



Figure [9]^I) presents the conductivity for different densities of states N = 10 9 - 10 11 
eV _1 m _1 as a function of the applied electric field. For low up to moderate electric fields 
the conductivity follows nicely the F 2 -behaviour, as we expected from the analytical 
expression previously reported [211 EE]. Specifically, when the condition F 2 /g(T) <C 1 
is satisfied, i.e. ea~ x F <C 2k B T, two of us have showed [26] 

hV' cr (F,T) oc hV' cr (0,T) + h(F)/f' cr (T), (37) 

where hW' cr (0,T) = -(T<j' cr /T) 1 / 2 , h(F) = F 2 and f> cr (T) = \(T l > cr /T) 1 ' 2 /2g(T)}- 1 . 
Increasing the electric field the conductivity becomes independent of the temperature 
and follows a l/F 1 / 2 -behaviour (Fig. [9](II)). The linear fit is of the form y = A + Bx 
with A = -15.664±0.018, B = -11159.837 ± 142.759 and R = -0.99959. In the region 
between the F 2 and l/F 1 / 2 -behaviour, In a increases almost linearly with F. 

The influence of both F and T on o is shown in Fig. [101 where we depict a/ao versus 

= eF/2ak B T, for (a) A^ = 10 9 eV"^" 1 and (c) N = 10 11 eV^^" 1 . F = (5 x 10 3 - 

1 x 10 8 ) Vm" 1 and T = (10 - 150) K. The influence of F on a depends on T, and it is 
greater at lower temperatures, while at higher temperatures the influence of F decreases 
significantly. Comparing Fig. [10] with the corresponding one for high temperatures (cf . 
Fig. [5]) we observe that the dependence of a/ao on (3 = eFj^aknT for low temperatures 
is very strong especially in contrast to the corresponding variation at high temperatures 
studied earlier on Subsection 13.11 

Fogler and Kelley [22], Raikh and Ruzin [25], and Ma et al. [23] refer also to 
the transition of a ID disordered electron system from the ohmic to the non-ohmic 
behaviour, and in this respect their results are consistent with our results, for both low 
and high temperatures. For strong enough electric fields, Fogler and Kelley [22J, Pollak 
and Riess [46], and Shklovskii [47], claim that in strong electric fields only the forward 
hops need to be considered, in contrast to Apsley and Hughes [40] who integrate over the 
entire space. Our methodology follows in that sense Apsley and Hughes [40] summing 
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Figure 9. (I) Left Panel, a versus F varying the density of states, i.e. for cases 
(a) iV = 10 9 eV^m -1 , (b) N = 10 10 eV-W 1 and (c) N = 10 11 eV^m" 1 . 
F = (5 x 10 3 - 1 x 10 s ) Vm- 1 and T = (10 - 150) K. (II) Right Panel, a versus F~ 1 / 2 
for case (c) N = 10 11 eV^m" 1 at T = 10 K and F = (4 x 10 7 - 1 x 10 8 ) Vm _1 . In 
both panels a -1 = 2 A. 




Figure 10. er/<7o versus (3 = eF/2aksT for cases (a) iVo = 10 9 eV 1 m 1 and (c) 
iV = 10 11 eV^m -1 . a' 1 = 2 A. F = (5 x 10 3 - 1 x 10 8 ) Vm- 1 and T = (10 - 150) 
K. The insets present enlargements of ct/ctq versus eF/2afcsT. 



over forward as well as backward hops in our ID polaron system. Bourbie et al. [48J 
study the F-dependence of the hopping conductivity in disordered 3D electron systems. 
They propose -among other mechanisms- that a decreases with increasing F, when F 
is strong enough to affect the tunneling probability. This is due to the influence of F on 
the number of percolation paths, in the sense that increasing F certain paths become 
disallowed. This is in analogy with our discussion about the destructive role of F at high 
F and low DOS. Bourbie et al. [29] . taking into account that U F affects the effective 
dimension of the transport path, reducing it in the high-F regime to ID" , showed that 
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when F = 10 8 Vm" 1 , T « (200 - 330) K, a' 1 = 2 A and JV = 10 12 eV^m" 1 , the 
conductivity decreases with increasing temperature. This has been attributed to the 
competition between thermal-assisted and field-assisted hopping. We have obtained a 
similar behaviour for a when F = 10 8 Vm -1 , at high temperatures T = (160-300) K, 
a^ 1 = 2 A and N Q = 10 13 eV _1 m _1 (cf. Fig. [3]). Bourbie et al. have also included 
different forms of the DOS and mention that these different DOS lead to very similar 
F-dependence of a. D. Bourbie [30] also used some different values for the extent of 
the 3D electronic wave function arriving at the result that greater extent of the carrier 
leads to higher conductivity in analogy with our results for ID polarons. However, we 
underline that all the above works [221 SQl SSI S7J SHI EH [30] refer to electrons while 
we study polarons. Moreover, in our work we have scrutinized the importance of the 
magnitude of the density of states and the spatial extent of the localized electronic wave 
function (for arbitrary electric fields up to the polaron dissociation limit and for any 
"reasonable" temperature). Finally, in ID systems the ionic or protonic transport might 
play a role in some cases [3SJ Efll EH [52]. However, in the present manuscript we do not 
investigate such possibilities. 

4. Conclusion 

We showed that the strength of the localization which determines the size of the 
formed polaron along with the density of states are two key factors for the conductivity 
and its dependence on the electric field and the temperature either at high or at low 
temperatures. These aspects of small polaron hopping have been nearly ignored in the 
past. 

To accomplish our task, we developed a novel theoretical approach inspired by the 
eminent work of Apsley and Hughes [40J in combination with the GMCM [131 El [El 
|24"1 [26] and references therein. In addition, we combined analytical work with numerical 
calculations. In the present model the expression which determines the conductivity 
(cf. Eq. [33| depends on both the density of states and the extent of the electronic 
wave function. We varied the DOS by few orders of magnitude near values which are 
relevant to common ID systems [271 1211 [291 El] and the extent of the electronic wave 
function from 1 to 5 A, i.e. for reasonable values for common organic molecules [311 [32] . 
Although in the present manuscript we used for simplicity a constant density of states, 
it is evident from Eq. [31] that one could also try an energy dependent DOS via the same 
approach. We examined cr(T, F) for temperatures from 10 up to 300 K and up to the 
electric field values where polarons dissociate («lx 10 8 Vm -1 ). 

We showed the the electric field plays both a constructive role by offering energy 
for the polaron hops and a destructive one, in the sense that the stronger it is the more 
it forces the polaron to jump opposite to the F direction prohibiting forward jumps to 
neighboring sites. The relative strength of these two roles depends on the DOS and 
localization regimes. 

Our present method confirms that either for high temperatures or for low 
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temperatures, higher density of states leads to higher conductivity. This is done in such 
a way that is smaller for higher densities of states. Conclusively, augmenting DOS 
by few orders of magnitude increases the conductivity by many orders of magnitude. 

For high temperatures, for the smaller densities of states % < i.e. the conductivity 
is larger for lower F. This is due to the competitive role of the directionality imposed 
by the electric field and the temperature. This directionality affects destructively o 
when only few sites are available for the polaron i.e. for small DOS. We noticed that 
according to Eq. [T^and Fig. Hfl), the electric field affects the range between two sites in 
the "hopping space" both for the absorption and the emission branch. On the contrary, 
this effect does not appear at low temperatures, because in the corresponding expression 
for the range between two sites in the "hopping space" at low temperatures (Eq. [25] and 
Fig. Efl)), the electric field affects only the finite area of the absorption branch. We also 
noticed that the electric field plays a constructive role, too, due to its energy offer to 
the polarons. For "medium" DOS the available sites are numerous enough so that the 
directionality of F hardly affects the conductivity. For higher DOS only the constructive 
energetic influence of the electric field appears. Now % > i.e. the conductivity is 
larger for higher F. Finally, we observed that the temperature has a greater effect on 
the conductivity, the smaller the density of states is. 

Our results confirmed that either for high or for low temperatures the behaviour 
In a oc T -1 / 2 previously obtained by two of us [26J, following a different theoretical 
treatment and taking into account the effect of correlations, holds for low up to 
moderate electric fields. Moreover, for low electric fields the conductivity follows the 
F 2 -behaviour [24"||26]. and increasing F the conductivity becomes independent of T and 
it follows a 1/ F 1 / 2 -behaviour while in the region between the F 2 and 1 / F 1 / 2 -behaviour, 
In cr increases almost linearly with F. 

We examined the deviation of conductivity from its ohmic value under the influence 
of both the external stimuli F and T (introducing = eF/2aksT). This was done either 
for high or for low temperatures, and for different DOS. We showed that <r(/3) depends 
strongly on the value of the DOS, and either decreasing or increasing <r{0) could be 
observed. We noticed that the variation of cr/cr versus /3 is generally very small in high 
temperatures compared to the corresponding variation at low temperatures. 

Finally, we studied the conductivity for different values of the spatial extent of the 
localized electronic wave function in the range a -1 = (1 - 5) A. Our results confirm that 
more localized polarons exhibit smaller conductivity. Particularly, five times increase of 
a^ 1 lead to two orders of magnitude greater conductivity. Moreover, we showed that 
4| > for any T when a^ 1 = 3, 4 and 5 A, while §| < when or x = 1 A. For the case 
a^ 1 = 2 A we observed an intermediate behaviour. 

In summary, we proved that the size of the polaron and the density of states are 
crucial factors for the behaviour of the conductivity and its dependence on the electric 
field and the temperature either at high or at low temperatures. 
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Appendix 



The integrals l\ and I2, relevant at high temperatures, are given below: 



/■i(2<„+^) 

'. = EL, Af(£?)[l -/(£?)] 

0,7T J S^i"^"" 006 ") 



Jt nn 2 J 2 » 
1 + 4 COS # 



cos OdEj. 



0,7T J 3(- B i-/ 35R »" COSe> ) 



(38) 



(39) 



As the integrals I± and J 2 diverge at j3 — — 2/ cos#, we change variables and integrate 
over R'. Hence: 
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The integrals ii, I2, h, h relevant at low temperatures are given below. 

n 



h = 



*C -e* + e; 



h = 



E / N(E*)[l-f(E*)]M l nn cos9dE*. 

0,7T J -°° 

Y. E . - N(^;)[i-/(^)]dE;. 



cos 6dE*. 



h = 



E 

0,7T 



cos* 



A(K*)[1-/(K*)]^. 



(41) 

(42) 
(43) 
(44) 
(45) 



As the integrals I\ and I 3 diverge at j3 — — 2/ cos6>, we change variables and integrate 
over i?'. Hence: 

h = E / ^(< n -^(l+/3cos^)+^))[l-/(< n - J R'(l+/3cos^)+E;)] J R'cos^ J R'.(46) 
/ 3 = E / iV(< n - J R , (l + /3cos^)+ J E*))[l-/(C -R'(l+(3cos9) + E*)]dR'.(47) 
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